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Abstract 

We study cold nuclear matter based on the holographic gauge theory, where baryons 
are introduced as the instantons in the probe D8/D8 branes according to the Sakai- 
Sugimoto model. Within a dilute gas approximation of instantons, we seek for the 
stable states via the variational method and fix the instanton size. We find the first 
order phase transition from the vacuum to the nuclear matter phase as we increase 
the chemical potential. At the critical chemical potential, we could see a jump of the 
baryon density from zero to a finite definite value. While the size of the baryon in the 
nuclear matter is rather small compared to the nucleus near the transition point, where 
the charge density is also small, it increases with the baryon density. Those behaviors 
obtained here are discussed by relating to the force between baryons. 
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1 Introduction 



It is difficult to study dense nuclear matter from the 4D non-perturbative approach 
like lattice gauge theory due to the sign problem, which appears when we introduce 
the chemical potential (See for example [TJ [2] ) . Then it is challenging to make clear the 
properties of the nuclear matter from the viewpoint of holographic gauge theory. In this 
approach, baryons can be introduced as solitons on the probe flavor branes [31 IH El E] • 
For the simplicity, we consider two-flavor case, Nf = 2. The soliton, which carries unit 
baryon number, is given as the BPST instanton solution of the SU(Nf) YM theory in 
the flat 5D space-time of the probe brane pE]. The instanton configuration has a scale 
parameter which could be identified with the baryon size. In the case of flat space, the 
size is free. On the other hand, the geometry of the world-volume of the embedded 
flavored brane is deformed in the fifth direction. As a result, the size of the instanton 
is fixed by taking into account of this deformation. Actually, in [J], this size of the 
baryon is determined by minimizing the action of the probe D8/D8, which is expressed 
as a sum of the Dirac-Born-Infeld (DBI) and the Chern-Simons (CS) terms. The CS 
term is also necessary to introduce the baryon chemical potential in this approach with 
the instantons, because the instantons, i.e. baryons, and the U(l) gauge field dual to 
the baryon number current are connected via the CS term. 

In this context, several approaches to study the nuclear system have been performed 
by introducing the chemical potential of baryons [3 El El HD1 [HI H21 H31 H3]- In |S], 
the instantons are, however, used to give a delta-function type source at the bottom 
of the D8 brane configuration only in the CS term. Furthermore, D4 branes are added 
to cancel the singularity generated by the delta-function source at this point, then the 
V-shaped D8/D8-brane configuration is obtained instead of the smooth U-shaped one. 
As a result, the authors of [8] have observed a gapless transition from the vacuum to 
the nuclear matter phase at zero temperature. However this point is not in favor of 
the expectation from other many kinds of theories (See for example [T5]). 

On the other hand, the authors of [10] have preserved the flavor gauge fields in the 
DBI action as well as in the CS term, then they could find the first order transition 
at finite baryon density. However, the configuration used for the flavored gauge fields 
is obscure, so that it is difficult to see the relation between it and the well known 
instanton. Furthermore, the profile of D8 brane is not solved in the form of U shape, 
which should depend on the gauge fields configurations in the same DBI action. 

Here we introduce the explicit form of instanton solution, which is exact in the flat 
4D space (our three-space and fifth one), by keeping its size parameter. Then it has 
a smooth configuration in the fifth coordinate direction instead of the delta function 
form. Furthermore, the flavored Yang- Mills fields are retained as instantons in the DBI 
action of the D8 probe brane in its leading order, namely up to the square of the field 
strength. The latter point is crucial to find a gap of the baryon density at the transition 
point as in the case of [10]. Differently from |10] , in our approach, the flavored YM 
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field is solved by determining the remained size parameter of the instanton. After that, 
the physical quantities like chemical potential and baryon number density are obtained 
to search for the phase transition. 

As for the profile of D8/D8 brane, it is restricted here to the antipodal solution. 
Then, we can use a simple D8 brane profile |4], which is obtained without any instantons 
as an antipodal U-shaped configuration. In general, the lowest energy solution for some 
finite baryon density n is not equivalent to the simple solution mentioned above, since 
it is not antipodal. However, the solution is always approximated by the simple one 
with a negligible correction for any n. This implies that the simple solution is useful 
at any n as the D8 /D8 antipodal configuration in our analysis. This is the reason why 
we restrict the profile to the antipodal solution. Another reason is for the simplicity of 
the present analysis. 

Through our analysis, we find a first order phase transition, which is expressed in 
n — plane, where n and /x denote the vacuum expectation value of charge density 
(times some dimensionful constant) and the chemical potential. At the transition 
point, the baryon density jumps from zero to a finite value, which corresponds to the 
transition from the vacuum to a nuclear matter phase. By adjusting the parameters, 
at the transition point, we observe the baryon mass as [is ~ 2.3GeV which is rather 
large compared to the realistic neucleon mass. 

Furthermore at this transition point, the size of the baryon is rather small, but it 
increases as sJJL at large /i. As a result, at about n — 4, the nucleons begins to overlap 
each other, then we arrive at the region where the dilute gas approximation is invalid 
at large [i. On the other hand, however, we could consider that this is interpreted as 
a clue to a new phase of quark matter. 

The outline of this paper is as follows. In the next section, the dilute instanton gas 
model is set up. Then, the embedded solution of the D8 profile and other fields are 
given in the section 3. In the section 4, energy density of the system is examined by 
using a simple profile solution to find phase transitions. In the final section, summary 
and discussions are given. 

2 Dilute instanton gas model 

Here baryon is identified with an instanton [4]. Then we approximate multi-baryon 
system by the dilute gas system of multi-instantons. The system is described by using 
D8 branes with the following action. 
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2.1 DBI action 

For stacked two branes, 

S Dp = -T P J d p+1 ee-*StrL, (2.1) 
where p = 7(8) for type IIB(A), and Str denotes the symmetric trace of flavor U(2) 



i > 



L = ^-det(/ + /i 

(fo)ab = (p MNd a X M 8bX N + B a b + F a bj Tq 



B a b — BMNd a X M d X N 
F ab = 2-Ka'Fab , 



ba 



(2.2) 

(2.3) 
(2.4) 
(2.5) 
(2.6) 



where a, b = ~ p, M, N = ~ 9, and To, T\ are the unit, Pauli's spin matrices. The 
Str part is expanded as 

StrL = Str | v /-det(/ ) (l - ^rx 2 + l(trx) 2 + ■")}, (2-7) 

x = fo'fu (2-8) 

where 'tr' denotes the trace of the coordinate index. The linear term vanishes for Str, 
then it is dropped. The coordinates are set as 

and we make the ansatz 



F a b 


= 2d[ a A b] , 


B a b = , 


A b 


= A b (z)6° b 






= Q(x m - 


a m ,p)e ijk r k 


l)iz 


= Q(x m - 


a m , P )r\ 



(2.9) 
(2.10) 
(2.11) 
(2.12) 



where p (a m ) denotes the instanton size (position), ei 2 3 2 = 1, i,j = 1,2,3 and m 
1, . . . , 4, where x 4 = z. 

Then the above series of fi, by retaining its lowest order, are obtained as 



StrL = 2 v /-det(/o)( v l + -Q 2 [(£ n ) +G 11 G ZZ \), (2.13) 
Gab = G MN d a X M d b X N . (2.14) 

Hereafter we neglect the higher order terms of the non-Abelian gauge fields. 
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2.2 Q as a Dilute Gas of Instantons 



In order to see the baryon spectrum, Q is given as an instanton solution in flat 4D 
space {x m } [3], 

2o 2 

Q = o. (2.15) 

((x m - a m ) 2 + p 2 ) 2 V ; 

However, it is not a solution of the equation of motion given by the D8 brane action. 

We use it as a trial function which is supposed to be a solution of the system. This 

is confirmed in our formulation such that the size parameter p of the given instanton 

configuration is determined to satisfy the variational principle of the energy density up 

to the leading order of the expansion given in (12.131) . 

According to the above strategy, we study multi-baryon state by replacing Q(z) by 

the multi-instanton form with dilute gas approximation. Then we rewrite Q as 

Q 2 = y (216) 

where the overlapping between the instantons are suppressed in obtaining Q 2 . So we 
find the energy density in the flat 4D space as a sum of each single instanton, 



POO 

2iVj / dzq{zf 
Jo 



2tt 2 

^-N It (2.17) 



where 



f ~ 2{z 2 + p 2 f/ 2 (2 ' 18) 
for the case of dilute gas, where interactions between instantons are neglected. Then 
the result is given by one instanton "mass" times their number Nj. 

As for the size parameter p, we determine it by minimizing the embedded D8 brane 
energy as in [3]. We notice that there is another method giving Q including its z 
dependence by solving the embedding equations of motion [10]. We, however, solve 
only for the size parameter here. 

2.3 CS term 

Supposing a reduced 5D action, we consider the following form of CS term for Nj = 2 

Scs = ^ emi "' m4 / dx»dz~A Tr(F mim2 F m3m ,) . (2.19) 

For the instanton configuration, we have 

"\N r O 2 
Scs = ^i dx^dzA -^—. (2.20) 



1 



This term is then included in our calculation as the coupling of A and instantons by 
supposing Q = Q(z) as in the DBI term. 



3 Solution of embedded D8/D8 brane 

Induced metric for D8 brane is given as [31 0] , 

A/ 2 /4 4 \ 

dsl = (-^{z^dxV + g(r)-k^%z)dz 2 + k l '\z)dQVj , (3.1) 



k(z) = l + z\ e~* = -^PHz)-^ , (3.2) 



where k = U 3 = U\ K + Ukkz 2 and we set 



Ukk = 1, R 3 = l = irg4N c , X = 3^^, M KK = ^, (3.3) 

where Z s and g s are the string length and the string coupling respectively. Moreover, 
we give 

g(T) = l + z 2 k 1 / 3 r' 2 , T ' = ^^. (3.4) 

oz 

Note that r should have a periodicity r ~ r + 5r in order to avoid a conical singularity 
at z — 0. In the case of (13. ip . 5r is denoted as 2tt [3]. 
Then we have 



S m = -T P J c/ 9 Ce-*StrL , 

= -f^ / d^fc^^StrV-det^ + FMjv), (3.5) 

g i M ) N dx M dx N = k 1/2 {z)7 1 ^dx^dx v + g{j)k~ hl \z)dz 2 , (3.6) 

Using the above formula, we find 

5 D8 = -2k J d A xdzk b '\z)^ g(T)k~ l l 3 (z) - E 2 Z (l + ^Q 2 (V 1 + fc 1/3 ^(r) _1 



where = F 2 o = <9 2 A (2:) and Q = ttiQ- Similarly we obtain 



(3.8) 

f) SlimimrKr tsto nnf qiti 

4Ao' 



3% 5 

Hereafter we use Q instead of Q for simplicity. 



o2 \3 J\r 

£0* = -^ I d*xdzA (z)Q 2 . (3.9) 
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The effective action to be solved in this case is given as 

S = S DS + S cs = -2kV 3 J dtdzL Q , (3.10) 



Lq = k 5 '^g{r)k-^-ElQ l -QA,nql 1 (3.11) 
Qi = 1 + §nqg + k^g(r)' 1 ) , (3.12) 



2 

where n is the density of instantons defined by 

n = ^. (3.13) 

Note here that we replaced the previous instanton configuration q 2 to the following 
new form 

2 9 Try 

% = H7( 2^/2 > ( 3 - 14 ) 

» {{z - z Q y + p 2 ) 

which is reduced to q 2 for zq = 0. In any case, we take the position of the inserted 
instantons at the bottom of the D8 brane since they are pulled to the horizon by the 
gravitational force. In other words, we put the instantons at their stable position. 
On the other hand, they are distributed at random points in the three space of our 
Minkowski space-time. 

3.1 Equations of motion for the system S = Sds + $cs 

We suppose the form ( I3.14p for the instanton contribution, where the parameter p is 
remained to be determined. It is given here by minimizing the action. 

The equations of motion for A Q (z) and r(z) are given as 

d z \k 5/6 . E * Qx\ = 6nq 2 , (3.15) 



gi^k-Vz _ E 2 

d z ffc 5/6 i " V Qi) -3nd z (e/ 2 ^g(T)k-V3-E*^)=0. (3.16) 



lg{r)k-^ - E 2 Z 

At first, we solve for Aq(z), whose equation of motion is given as 

fc 5 / 6 . 5 ( 1 + Inq 2 (k~ l + k^/gir))) =d, (3.17) 



g ( T )k-i/z - El V 2 



d = 6n dzq + C = -nn (z — 20)77 + ^ ■ (3.18) 
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Here we take the integration constant C as zero since the displacement would vanish 
when instantons disappears (or n = 0). Then we have 

„1/2l-1/6 L5/6 

In the next, by substituting the above E z into (13.1 6ft . we find 

*Iv^ e , / ^L,, (3.20, 



9 1/2 M I / " vT + P 9(t) 

where Co denotes an integration constant, which is related to the bottom coordinate z . 
The solution r{z) has in general has a U-shape configuration in the r-z plane. And, 
we find t'(zq) = oo at its bottom point, then the following relation is obtained from 

z Qi(z )k 3/A (z ) = c . (3.21) 
3.2 The simple solution in the antipodal case 

Here we consider the antipodal solution which satisfies r(oo) = tt/2 (—it/2) in the limit 
zq — > 0. In this case, we obtain a simple solution of Eq. (l3.20p as 

t\z) = for z > zq. (3.22) 

This is the exact U-shaped solution at n — 0. 

In general, the solution of Eq. (13.201) depends on the instanton density n, and we can 
see it by solving equation expanded perturbatively by the power series of n. In the 
case of zq — > 0, we can obtain a solution^ which behaves r(oo) > tt/2. However, such 
D8(D8)-branes configuration at zq — > degenerates and can have any value of r(oo), so 
we can adjust any degenerated solution to be the antipodal solution (I3.22p by adding 
the part of the solution which does not contribute to the action (I3.10p at z — > 0. In 
this sense, the value of r(oo) is not just tt/2, but tt/2 + 2nm. 

On the other hand, we can consider other antipodal solutions that satisfy r(oo) = 
7r/2 at finite zq and n. However, it is confirmed by numerical analyses that the energy 
of the solution at finite z is always larger than the solution (I3.22p in the same n. 
Therefore, we adopt here the solution (I3.22p . 

Then E z and the chemical potential are obtained by noticing g(r) — 1 as 

E z = 3.23 

>Qlk 5 / 3 + d 2 



1 It seems to apparently have a conical singularity in the limit zq — > when we see it in polar 
coordinates (t,z). However, the solution has no singularity even if we take the limit since we impose 
the boundary condition t'(zq) = oo for the solution to be smooth. 
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Qo = l + lnq^k^ + k 1 / 3 ), (3.24) 



u = n c + dzE z . (3.25) 
Jo 

Here the quark chemical potential \i and n c are defined as 

fi = A (oo) , // c = 4,(0) . (3.26) 



E* = — (-,--, + ■■■) , (3-27) 



In this case, -E 2 is expanded at large z as follows 

9mr 2 / 1 1 

Here we notice the coefficients of the leading and sub-leading orders are independent 
of p. Then we obtain 

9nir 2 1 

Ao = fi-— - + ••• . (3.28) 

The charge density is expected to correspond to the coefficient of the sub-leading term 

as 

c^W) = n = -y- , (3.29) 

where c is a dimensionful coefficient obtained by returning dimensionless variables to 
dimensionful ones. 

In the small z region, E z is expanded as 

„ 9n7r 2 / o \ , 

£ 2 = — (a x z + a 2 2 3 + •••), (3.30) 

where ai^ are the functions of p. Then we find the asymptotic solution of A Q , 

9nir 2 a-i 9 

A = fi c + ^^z 2 + -- - . (3.31) 

where /i c is undetermined. 



4 Energy density and phase transition 

The Chern-Simons term is rewritten as 

/*oo /*oo 

- 6n / dzAog 2 = -6n / dzA (d z q b ) (4.1) 

JO JO 



/ f°° 07T \ 

6n(/ dzE g q b -—fi\ (4.2) 
6n f jH (q b ~^f\- , (4.3) 



2 We notice here that the chemical potential given above expresses the one of the quark rather than 
that of the baryon. 
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where we used qb(oo) = vr 2 /3, g&(0) = and Eg. (13.251) . We notice about the sign 
of the Scs term. When it is taken as minus as given above, and then we use the 
self-dual instanton solution. In this case we find positive chemical potential as shown 
in Eq. (l3.25p . On the other hand, for opposite sign and self-dual solution, we obtain 
opposite sign chemical potential 

POO 

fi = ii c - \ dzE z . (4.4) 

JO 

and the same sign energy density except for the chemical potential part, 

r°° ( r°° ( tt 2 \ tt 2 \ 

+ 6n J dzA q 2 = 6n I J dzE z I q b - — J + — p c J . (4.5) 

When we use anti-self dual solution, then the same expression is obtained for the 
opposite sign. So we must use the sign in our case, we should use the self-dual solution 
(baryon). When we apply the same analysis for the anti-baryon, we must change the 
sign of Scs and use the anti-self dual solution. Then we find the same result in both 
cases. 



By using the above solution for Aq, we define the energy density as 

E{n) = V 3 (I dzL Q - nfi c ^j , (4.6) 

where Lq is given as 

Lq = Ldbi + Lcs , (4.7) 

Ldbi = , ^ , (4.8) 

Ql + ^-5/3^2 



37T 2 



'CS 



6nd . Qb 4 , (4.9) 



kjQo + k-V 3 d 2 



3ir 2 2z 3 + 3zp 2 

» - (4 ' 10) 

Then, we compare the energy density E(n) with E(0) to find the critical value of 
chemical potential, \i cr . Above which we will find E(n) < E(0) then the baryon phase 
is realized there. 

The integration of the Lagrangian Ldbi for finite n diverges, so we need a subtrac- 
tion, which can be chosen as the energy density of the vacuum with n — 0. E(0) is 
given as 

E(0) = V 3 Jdzk 2/3 . (4.11) 
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Then we can perform the calculation of E(n) — E(0) by using Ldbi, 

Ldbi = L DBI — k 2 ^ 3 . 



We can thus estimate the energy difference for a given density n by normalizing as 



In estimating this energy density, we set fi c = A o (0) = since A (z) must be a two 
valued function which is symmetric at z — 0. Performing the integration over z for a 
fixed n, we obtain S(p) as a function of p. 

In the estimation of S(p), we must be careful about the value of £(0). The value 
is sensitive to the order of taking the limit of z — > zq and p — > 0. Our estimation is 
performed by the integration with respect to z from Zq to z = oo, then the limit of 
p — > is examined. 



Fig. 1: £{p) — £(0) versus instanton size p near phase transition point p cr ~ 1. 



4.1 Phase Transition 

To estimate physical quantity, we must determine the size parameter p by minimizing 
the energy density S(p). This is performed by estimating S(p) for each fixed p, not 
for fixed n since the parameter which characterizes the boundary theory is not n but 
the chemical potential p,. In our formalism, however, the energy density is given as a 
function of n and p from the Eq. ( 14. 121) . At the same time, p, is also obtained as a 
function of n and p from Eq. (13. 25ft . Then 8{p,p) for a fixed p, is estimated by using 
both calculations of S(p,n) and /x(p, n) although it is a little complicated. 

Noticing the property of £{p) mentioned above, we estimated £{p) for a fixed /i. 
Then, fix the the value of p as p m at the minimum of energy density, namely as 




(4.12) 




O.d 
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£ m in = £(Pm) f° r the given p. After that, we determine n. The resultant curves of 
£{p) are shown in the Fig. [T]near the phase transition point. The phase transition is 
observed at p cr = 1.0176, where £(p m ) changes its sign from positive to negative. This 
implies the transition from the vacuum to the nuclear matter phase at a finite charge 
density and finite chemical potential, (p,n) = (1.01,0.85). This transition is also seen 
more clearly from the Fig. [2] represented in the plane of p-n. 

The reason why n is finite at the transition point is reduced to the preserved flavored 
gauge fields in the DBI action, this is reflected to the factor Q > 1 given in (]3.24p . This 
point is also mentioned in the introduction. In our model, this is explicitly understood 
as follows. When we set Qq = 1, the integrand in (14. 12p is expanded near n = as, 

Ldbi + L cs = 6nd qo ~ 3 k ni/A - ±-d 2 + ■ ■ ■ , (4.13) 

where • • • denotes the higher order terms of n, and notice d is the order of n. It is easy 
to see that the first term of the right hand side of (I4.13P is non-positive since 

3ir 3 

qo- — <0. (4.14) 

Here the equality is obtained for p = 0. Then the energy density is always negative for 
n > in this case, namely for Qq = 1. This implies the transition occurs at n = when 
we neglect flavored gauge fields in the DBI action. So we can say that the interaction 
among the baryons as dilute gas is included in the model through the preservation 
of the flavor gauge field configuration in the probe brane. In spite of the dilute gas 
approximation adopted, the attractive force to form the nuclear matter phase would 
be reduced to the curved space in the fifth coordinate direction, namely the interaction 
reduced to the vacuum of the Yang-Mills theory. Although the details of the dynamics 
should be examined more, the importance of the flavored gauge fields in DBI is clearly 
understood from our model setting. 

Size of baryon 

As for the size of the baryon, which is identified with p m , its p dependence is shown 
in the left of the Fig. [3j From this numerical results, we can see that p m grows with p 
monotonically, and it is fitted by the curve 

p rn ~ yjfj,- 0.78 -0.10. 

This implies p m increases with the increasing baryon density since p grows with n. This 
behavior could be naturally understood by considering the attractive force between the 
baryons to form the stable nuclear matter. This force becomes large with increasing n 
since the number of the nearest baryons of one baryon grows. As a result, the attractive 
force between them increases, and then the size of the baryon is stretched to larger 
size. 
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However, at some critical density, baryons begin to overlap each other, then it would 
become difficult to discriminate each baryon spatially. In this case, we will find a new 
nuclear phase, but our model cannot be applied to this density since the dilute gas 
approximation is no longer useful in this region. In the right of the Fig. [HI this point is 
seen by considering a cube surrounding instanton in three dimensional space which the 
length is equal to 2p. We can see from this figure that the dilute gas approximation 
is not valid any more for \i > /xi imit . We notice here the following point related to the 
parameter A. Since the critical value of \i depends on A when we adjust the parameters, 
we find that the point /Unmit increases as A deceases. In the case of A > 155, our model 
would be no longer valid since p\\ m \t become smaller than fi CI . 




Fig. 2: The diagram of (/i, n = ^f-ra) for the simple solution. The first order phase transi- 
tions, from vacuum to the nuclear matter phase at fi cr = 1.0176, is seen. 



Numerical estimation of the size 

In the above analysis, the physical quantities are set to be dimensionless. In this 
setting, we get 



Her ~ 1.01 , n cr ~ 0.019 (n cr ~ 0.85) , p m ~ 0.38 



(4.15) 



at the critical point. So we must see their dimensionful values to compare them with 
the realistic physical values, (/ie = N c /i, n B , p B ). 
In our notation, furthermore, we have set a^H 



1 27tt A A . A2 /9tt\ 2 
2N f 2A V A / 



(4.16) 



3 In reality, we dropped the factor 2ira' in (|2.15[) . However this factor can be absorbed into n 
because they always arise together. So we can recover this factor by interpreting n as (27ra') 2 n. 
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Pm 




Fig. 3: (Left) The Plot of (//, p m ). The solid line is given as p m = y/p - 0.78 - 0.10. 
(Right) The relation p-p s n for A = 16.6. The horizontal line at p 3 re = 1/8 corresponds to a 
rough accommodation limit of the dilute gas approximation. The intersection point at about 
fi ~ 3.1 denotes the useful limit of the model. 



Therefore when we compare p cr , n cr , and p m given above with physical quantities of 
baryon, their scales are given by Mkk and A from the following relations as 

WfN&r^- ~ 6 x 1.01 x (4.17) 



M KK v * * 27tt 27tt 

1 1 / A \ 2 

J^L- = n -~ 0.019 x 7 - = 0.019 — , (4.18) 

m kk (lira 1 ) 2 {2na') 2 \ 97r / 

PbM kk = p m ~ 0.38. (4.19) 

Here we used Nf = 2 and N c = 3. 

In the Sakai-Sugimoto model[5J [S], these parameters are chosen as M^k — 0.949 
GeV and A = 16.6 to fit the values of p meson mass and the pion decay constant in 
the vacuum. Then these parameters lead to 

He ^ 2.3 GeV, n B ~ 0.70 fm" 3 , p B ~ 0.079 fm. (4.20) 

This value of % is close to 4no where n = 0.16/m -3 is the density of the normal 
nuclear matter. It is interesting that 4uq corresponds to the density of the inner 
core(n > 2no) of a neutron star [TBI . 



5 Summary and Discussions 

We have studied cold nuclear matter using the Sakai-Sugimoto model, where the baryon 
is introduced as a soliton on the probe flavor branes. The soliton, which carries unit 
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baryon number, is given by the BPST instanton solution in the probe D8 branes. Multi- 
baryon system is approximated by the dilute gas of instantons in the brane action. The 
profile of the D8 brane depends on the configuration of the instanton and its density. 
So it is rather complicated to solve them from the given D8 brane action. 

The instanton configuration used here has the same form with the one given in the flat 
4D space, but it is not the solution in the deformed bulk background. So we recovered 
this point by adjusting the size parameter of the instantons, which is identified with 
the baryon size. As for the D8 brane profile, we solved it numerically and found that 
the lowest energy configuration is well approximated by a simple solution, which is 
obtained without instantons for the antipodal configuration. Through our analysis, we 
have restricted to this configuration as a boundary condition of the D8 brane which 
fixes the dual theory. 

In our model, we could find the first order phase transition from vacuum to nuclear 
matter accompanied by a certain jump of the baryon density. This jump could be 
interpreted as the effect of the attractive force acting among the baryons since a finite 
density of the instanton gas in the DBI action is needed to form a baryonic system. 
Actually, we cannot observe the jump in the models which lack such a term in the DBI 
action. While the density is rather large compared to the nucleus, our numerical results 
obtained near the critical point seems to be reasonable. As a bonus, we could find the 
size of baryon for each /x. It is rather smaller than the nucleon size just after the 
nuclear matter is made, but it increases with the density or chemical potential. This 
observation would be reduced to the increasing attractive force among them. Then we 
expect a new transition from the nuclear matter to the quark matter at enough large 
density, but it is out of our dilute gas approximation. 

We should notice on the freedom of parameter choice in determining the physical 
quantities. For example, the numerical values of n B and pB obtained here depend on 
our parameter choice. If we set Mkk — 0.496 GeV and A = 61.7 [32], then we get 

/ie~4.4GeV, n B ~1.4fm" 3 , p B — 0.15 fin. 

This value of n B is much denser than the ordinary nuclear matter. In this case, we 
obtain 8p B ns = 0.058, so the nucleons don't still overlap each other in spite of the 
high density. Then it is still useful to apply the dilute gas approximation in this 
parameter choice. Therefore it would be possible to extend our analysis given here to 
high density nuclear matter like neutron stars. Another impotant extension is to the 
finite temperature YM theory by retaining the chemical potential. Then we could see 
richer phases of nuclear matter. This point will be visited in the near future. 
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